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Dalitz Analyses in Charm 
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Dalitz analyses in charm give access to magnitudes and phases of the charm decay amplitudes. They play a 
significant role in charm-mixing measurements, in the measurement of the CP-violating phase 7 in B decays to 
charm, and in the analysis of light meson resonances. We review recent results in all three categories. 
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1 . Dalitz Analyses 

Dalitz analyses give access to the full complex decay 
amplitudes, allowing the measurement of magnitudes 
and phases. This can be exploited in several ways: 

• Explore the properties of light meson resonances 
such as the p, ai, f(980) and more controversial 
ones such as a and n. 

• Investigate charm itself, in particular mixing 
and CP violation in charm. 

• B mesons decay to charm most of the time. The 
phases of the B decay amplitude can in cer- 
tain circumstances be measured by analysing 
the charm Dalitz plot. This provides the the- 
oretically cleanest and statistically most power- 
ful direct constraint on the CKM phase 7 that 
is currently known. 

In the following we will consider each of these items 
in turn. 



with A = (rn\) — m\ — m^)'^ — Arn^rn^g within the 
kinematically allowed limits, and A = outside. In 
the above expression, OiC*^" describe complex contri- 
butions to the total decay amplitude. In the simplest 
case, aie"*' are complex Breit-Wigner distributions (or 
similar e.g. the Flatte distribution 0]) describing indi- 
vidual particle resonances, with additional factors tak- 
ing into account angular momentum conservation, and 
form factors (Blatt-Weisskopf penetration factors Q). 
This so-called isobar model has some shortcomings, 
the most severe one being that it violates unitarity, 
especially in the case of wide, overlapping resonances. 
More complicated models such as the K-matrix for- 
malism [11,1^, which respects unitarity, may therefore 
be necessary to adequately describe the observed data, 
and to provide a theoretically satisfactory model. The 
general consensus - at least amongst experimentalists 
- appears to be that the isobar description is adequate 
for P and D wave resonances, but not for wide S wave 
resonances. The adequate description of i = decays 
is one of the main topics of interest in the next section. 
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2. Dalitz Plots 

The kinematics of a 3 body decay D ^ A, B, C 
(such as D+ K+K^7r+) can be fully described by 2 
parameters. In terms of the four momenta of the three 
decay products, which we will denote as pa,Pb,PC', 
one usually picks the following invariant-mass-squared 
parameters: 
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These parameters are Lorentz invariant, and the 
phase space density in terms of these parameters, 
—7 — „ \ f, — 5 — r is flat, i.e. constant inside the kine- 

matically allowed limits, and zero outside. A Dalitz 
plot is the decay rate in terms of these or equiva- 
lent variables, displayed in a 2-dimensional plot. The 
full decay rate is given by [l| : 

d^T 
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3. Dalitz Analyses and Light Meson 
Resonance 

3.1. The K, a problem 

The S-wave resonances a tt+tt^ and k K+tt^ 
are needed to describe the data in isobar fits to 
D+ 77+77-7:+, DO Ks7r~7r+ and D+ K-7r+7r+. 
However, it is unclear if this is compatible with 
LASS scattering data, and fits to D° K"7r+7rO, 
D^ 7r"*'7r~7rO do not require the addition of a or k. 
K-matrix models do not explicitly require cr or k. 

A wealth of measurements and interesting informa- 
tion has been published on this topic. Here, we will 
only consider two decay channels for which we have 
recent results, one for the tt+7t~ S wave (a), and one 
for the K+7T S wave (k). 
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Recent analyses of this channel include E791's anal- 
ysis using an isobar fit with a a resonance , and 
FOCUS, who pioneered the K-matrix approach in 
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Figure 1: The plots show the magnitude and phase of the 
amplitude of the S wave component for the different mod- 
els considered in CLEO-c's D° — > tt+tt"?!"'" analysis 
The dotted lines represent the upper and lower limit of 
the range of values obtained for the various isobar mod- 
els. The broken line corresponds to Achasov's model, the 
solid line to Schechter's model. There is good agreement 
between the models. 



3.4. Four-body "Dalitz" Analysis 

Essentially the same formalism as for 3 body de- 
cays can be applied to to 4 body decay. Such anal- 
yses are challenging as the equivalent of the Dalitz 
plot now has 5 dimensions instead of 2, phase space 
is not flat in the usual invariant-mass squared vari- 
ables, and the amplitude structure is more complex. 
A recent example of such an analysis, using an iso- 
bar model, is given by FOCUS for the decay chan- 
nel D" TT+TT+TT-TT- FOCUS obscrvc that 
ai(1260)7r is the dominant decay channel, fol- 
lowed by D° pp. The authors find that the ai 
predominantly decays to an. Many more results can 
be found in the paper, including the pp polarisation. 
Four body amplitude analyses of D° decays also play 
as significant role in extracting 7 from — > D°K='= 
decays, which is discussed in section [5] 



4. Charm Mixing and CP violation with 
Dalitz Plots 



this channel, and also analyse Dg ^ 7r+7r^7r+ in their 
study 0]. The most recent result is by CLEO-c 
using ~ 2600 signal events. CLEO consider isobar 
models with different descriptions of the fp (980) and 
(7, and two models that respect unitarity and chiral- 
ity, one according to Schechter and another de- 
veloped by Achasov [ll| . All models considered agree 
with each other and the results are consistent with 
previous fits. The amplitude and phase of the S wave 
contribution as function of tt+tt^ invariant mass is re- 
produced in figure Fig[TJ 



3.3. D 
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The branching fraction of D+ — > K 7r+7r+ is com- 
parably large, with (9.51 ± 0.34) % Over 60% of 
its decay rate proceeds via a Ktt S'-wave, as has been 
observed in several experiments. In 2002 E791 [isj . 
using an isobar model, found a large k contribution. 
In 2006, E791 re-analysed the same data with a model- 
independent description of the S wave, using a binned 
amphtude and phase [3. In 2007, FOCUS [3 ap- 
plied the K-matrix formalism, constrained by LASS 
scattering data to 54fc events. The most re- 

cent result, which also has the largest data set, is 
from CLEO-c in 2008, using 140k events, with very 
little background (1.1%) CLEO-c fit the data 

using both the isobar and the model-independent ap- 
proach, and compare their result to models used by 
other experiments. For both types of model, CLEO- 
c get a significantly improved fit if they allow for 
an isospin=2 tt+tt" S wave contribution, where the 
model-independent approach gives the better pcr 
degree of freedom. 



4.1 . The neutral D° system 

The neutral D system is the only neutral meson 
system that mixes and consists of up-type quarks. 
It therefore provides a unique window on Flavour 
Changing Neutral Currents (FCNC's) between up- 
type quarks, which can be affected by new physics in 
a very different way than those in down- type quarks, 
such as investigated in the ^ systems. Particularly 
interesting is CP violation in the D system, which is 
essentially zero in the Standard Model, but could be 
significant in many New Physics scenarios. The dis- 
covery of CP violation in charm decays would be a 
clear signal for New Physics. 

A. Mixing Parameters 

The physical mass eigenstates Di, D2 of the neu- 
tral D meson arc superpositions of the flavour-specific 
states D° and D": 



Di = pT)° + qT)° 
D2 = pD°-gD° 



(4) 
(5) 



where p and q are complex numbers satisfying \p\'^ + 
Igp = 1. The mass and the width difference between 
Di and D2 are Am and AF. Mixing in the neutral D 
system is conventionally parametrised by the param- 
eters X and y given by 
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In the absence of mixing, both parameters are zero. 
The Standard model expectations for x and y are 
10-3 - 10-2 [H. 



B. CP Violation in Charm 



If 



7^ 1, CP is violated (CP violation in mixing). 

CP violation in the interference between mixing and 
decay is parametrised by the phase 0. The phase is 
the exact equivalent in the neutral system of the 
parameter —2/3 in the neutral system measured 
in B*^ — > J/t/iKs by the B factories. The phase of the 
ratio ^ is convention dependent and without fixing the 
convention, its value does not say anything about CP 
violation. Usually, however, a convention is chosen 
where 3. = 1 q^4> gg that both kinds of CP violation 
discussed here are encoded in the same complex ratio 

q/p- 

In the SM, CP violation in charm is zero for all prac- 
tical purposes, i.e. relative to the sensitivities of cur- 
rent experiments and those planned for the foreseeable 
future. However, many NP models predict CP viola- 
tion in charm at a level that would be experimentally 
accessible. A review with further details on charm 
mixing, CP violation and its sensitivity to physics be- 
yond the Standard Model can be found in [l^ . 



BELLE collaboration analysed approximately 0.5 M 
D°,DO ^ KsTT+TT- events, with the result 

X = (G.81±0.2Gt°:}?)% 

y - (0.37 ±0.251°:}°)% 

\p/q\ - 0.86±0.30t°:i° 

6 = -14° ±18° ±5° 



where the first error is statistical and the second error 
is systematic. This measurement dominates the world 
precision on x [1^. 

4.2.2. CPV measurements in SCS decays 

The statistically dominant measurement of CP vi- 
olation in charm from D° Ktt and the Dalitz plot 
method with D° TTTT study CP violation in ampli- 
tudes where doubly Cabibbo suppressed amplitudes 
interfere with D° mixing and a subsequent Cabibbo 
favoured decay. However, decays involving singly 
Cabibbo suppressed decays could in principle be af- 
fected by New Physics in a different way and hence 
show CP violation even though CF and DCS modes 
do not. 



4.2. D mixing using Dalitz Plots 

Evidence for charm mixing has been observed first 
using the "wrong-sign" decay D° K+tt- [2l|, 
[2^ . The final decay rate has contributions from the 
DCS amplitude D° K+tt" and a combination of 
mixing and a CF amplitude, D° — > D° ^ K~^tt~. A 
time-dependent measurement is sensitive both to CP 
violation and to the mixing parameters x'^ and y' . 
The primed parameters are related to x, y: 

[ ^' \ = [ cosJa'tt sin^A'TT W ^ ^ (g) 
\y' J \- sin Sktt cos Sk-^ ) KV ) 

where Sn-n is the phase difference between the am- 
plitude D° K+TT^ and D° K+tt-, measured at 
CLEO-c to be cos(5k^ = 0.9 ±0.3 [H. 

4.2.1 . D mixing and CPV in D ^ KsTttt 

CLEO-c pioneered the Dalitz plot analysis of the 
self-conjugate decay D° — > KsTt+tt^ for the D° mix- 
ing analysis [23|. In this case, CF modes such as 
D0^K*-7r+ and DCF modes like D° ^ K*+7r-, 
which contribute to the mixing measurement in a simi- 
lar way as in the two-body case, are in the same Dalitz 
plot, so their relative phase can be measured. The 
method gives direct access to x and y, and is also sen- 
sitive to the CP violation parameters \p/q\ and (f). The 



A. D" ^ TT+TT-Tr" and D" ^ K+K'ttO 

BaBar measure several CP violation observ- 
able the 3-body d ecay s D°/D° tt+tt-tt^ and 
D°/D° K+K-7r° df. BaBar perform a 

model independent binned analysis of the rate 
of D° 7r+7r-7r° vs D° ■k^'k~'k'^ across the Dalitz 
plot, as well as a model-dependent one, comparing 
magnitudes and phases of the CP-conjugate ampli- 
tudes. No significant CP asymmetry has been found. 
The authors conclude that any CP violation in the 
singly Cabibbo-suppressed charm decays occurs at a 
rate which is not larger than a few percent. 

BELLE consider the asymmetry of the total decay 
rate in D°/D° 7r"'"7r-7r°, integrated across the whole 
Dalitz space). Their result [111, Acp = (0.43 + / - 
1.30)%, also shows no evidence for CP violation in 
SCS decays. 



B. D+ ^ K+K-7r+ 

CLEO-c fit an isobar amplitude model to the Dalitz 
plot of the SCS decay D+ K+K-7r+ and its CP con- 
jugate, D- — !■ K-K"''7r-, and form the asymmetry of 
the decay fractions (oc |Ap) for each of the ampli- 
tude contributions. No evidence for CP violation is 
observed. 
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5. Dalitz Analysis in Cliarm for Precision 
B physics 

5.1. Introduction 

A central aim of current and future flavour physics 
experiments is the precision determination of the CP- 
violating phase 7. In terms of the elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) quark mixing 
matrix, 7 is defined as &-i:g{-V*fyud/V*^Vcd)- 

A theoretically clean and statistically powerful 
method exploits the interference of D°K* and 

— > D°K^ decays, where the D" and tf decay to 
a common final state / [1^, [10, Hf], [H, HI]- Suit- 
able final states / include 2-body states such as KK, 
TTTT [2^, [sO], Ktt [3^. 3-bodv final states such as 
KsTT+TT" and Kttt t" l32l . [33| and 4-body final states 
such as KirTTTT [3l|, |34| and KKtttt [11] . 

All such measurements are sensitive to the ampli- 
tude ratios 



rse 



i(5B±7) 



(9) 



In all cases, the measurement is affected the properties 
(especially phases) of the D° decay amplitudes. This 
is where charm physics can make a significant con- 
tribution to precision B physics. By measuring the 
phases of the charm decay amplitudes, the uncertain- 
ties in the 7 extraction in B physics experiments can 
be significantly reduced. 

5.2. CLEO-C and B± ^ D(Ks7r7r)K± 

The best direct constraints on 7 come from mea- 
surements in B^ ^ D{K^'k)K^ and related modes, 
at the B factories [s^ [lo]. The combined result is 
7 = 67°l25° El- dominant systematic uncer- 

tainty in these measurements is the model uncertain- 
ties in the description of the decay amplitude, cur- 
rently between 5° and 9°, which would soon limit the 
precision of this measurement at the next-generation 
flavour physics experiment, LHCb. 

CLEO-c's quantum correlated DD pairs give model- 
independent access to both magnitude and phase in- 
formation of the decay amplitude across the Dalitz 
plot. This additional information can be used as in- 
put for a model-independent extraction of 7 from a 
binned Dalitz plot analysis [13, [13, [H, [13] , thus elim- 
inating the model-uncertainty. 

As Dalitz plot variables, we use the invariant-mass 
squared of the Kg7r~ and the Kstt^ pairs, which we 
denote as s_ and s+ respectively. The phase of the 
D decay amplitude at a given point in Dalitz space is 
6^'=^'^^ {s-, s+). For the phase-difference between the 
D KsTTTT amplitude and the D — ^ Kstttt amplitude 
at the same point in Dalitz space, we define 

A,(s_,s+) ^ <5^---(s_,s+) - (10) 




0.5 



1.S 2 2.5 
m2(Ksn*)/GeV2 



Figure 2: Binning used for this preliminary CLEO-c re- 
sult. The binning is uniform in lAa], with bin 1 centred 
at 0°. 



The quantities measured by CLEO-c that provide the 
input to the 7 analyses of the B-factories and LHCb, 
are the averages of cos A^^'^'^ and sinA^^'^'^ for each 
bin, Cj and s,: 



a = (cosAi-)^ 
Si = (sinAa),; 



(11) 
(12) 



where the index i denotes the ith bin. The analysis of 
Kltttt events in the similar way to Kstttt, provides fur- 
ther input to the c,; and Si measurement in the KgTTTr 
Dalitz plot. The clean environment at CLEO-c al- 
lows the Kl reconstruction from kinematic constraints 
with high purity. The choice of binning affects the 
statistical precision of the analysis - it is beneficial to 
choose the bins such that the phase difference A^ , de- 
fined in Eqlini varies as little as possible across each 
bin (37! . The binning used for the preliminary CLEO- 
c results presented here is based on the BaBar isobar 
model [38| . A uniform binning in As , with eight pairs 
of bins (arranged symmetrically with respect to the 
diagonal axis defined by s_ = s+) is chosen. This 
binning is shown in Fig [21 The latest preliminary 
CLEO results for and Si from the combined anal- 
ysis in both KsTTTT and Kltttt Dalitz plots are shown 
in Table [H When used as input to the 7 extraction 
in the Kstttt mode at the B factories and LHCb, this 
is expected to replace the current model uncertainty 
of 7° — 9° with an uncertainty due to the statistically 
dominated error on c, and s, of about 1° — 2° [431. 



5.3. Coherence Factor 

The decay rates B^ D(hh')K*, where hh' stands 
for any two-body final state accessible to both D° and 
D°, are sensitive to 7 [l^, [13, [HI. For example for 
D" ^ K+TT": 
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bin number 


Ci 


Si 


1 


0.742±0.041±0.022 


-0.022±0.168±0.096 


2 


0.607±0.073±0.038 


0.009±0.220±0.054 


3 


0.064±0.077±0.052 


0.548±0.198±0.096 


4 


-0.492±0.133±0.056 


0.124±0.227±0.074 


5 


-0.918±0.053±0.039 


-0.118±0.194±0.057 


6 


-0.743±0.071±0.033 


-0.296±0.203±0.070 


7 


0.156±0.092±0.050 


-0.870±0.183±0.062 


8 


0.398±0.047±0.020 


-0.438±0.146±0.041 



Table I Preliminary CLEO-c results for the measurement 
of Ci and Si in the D Kstttt Dalitz plot, using input 
from both D Kstttt events and D — > Kltttt events in 
818pb"\ The first error is the combination of the sta- 
tistical error and the uncertainty that arises from the use 
of Kltttt results for the Kstttt Ci and Si determination. The 
2nd error is the remaining systematic uncertainty. 



T{B-^ {K+tt-)dK~) cx 

rl + (r-f + 2rBrf - cos (,5^ + 5^^ - l) , (13) 
where we used, in analogy to Eq[5] 



A{W K+Tf- 



_ Kit MS^'') 



(14) 



The addition of 3 and 4-body decay modes of the 
D such as D'' K+tt^tt" and D° K+tt^tt+tt^ can 
significantly improve this measurement. In multibody 
decays D° — > / , the resonant substructure needs to be 
taken into account, which can be achieved by adding 
only one additional parameter to describe the decay 
rate [U- For a generic final state /: 



T{B~ ^ {DdK-) (X 

2rBRfAfAf cos (Sb + 5^d - l) (15) 



where Rf is the coherence factor which satisfies 
< i?/ < 1; is the average strong phase difference. 
The larger Rf, the higher the sensitivity to 7 in a 
given mode. Decays of quantum-correlated DD pairs 
at CLEO-c can be used to measure both Rf and 8{j. 
CLEO-c measure the following values for the coher- 
ence factors in the K+tt — tt+tt" and K+7r~7r° mode, 
as presented for the first time at CKM 2008 



Rkstt 



On 



< 0.57 at 95% CL 

0+85° 



R 



161' 



48° 



0.79 ±0.08 



197°i?„*° 



The confidence regions in i?^ and 6^ space for f = K37r 
and f = KttttO, as presented at CKM 2008 are 
shown in Fig [31 The small coherence factor in K37r 
implies that this mode on its own would not be 






X' 'DOF =0-9/\ 




CLEO-c preliminary B16 pb 





Figure 3: 1,2 and 3 a confidence regions in Rf — Sj^ space 



for D' 







K^ 



V ■ 

and K+tt"tt". The star rep- 



resents the central value. Preliminary CLEO-c result, first 
presented at CKM 2008. 



very sensitive to 7. However, this coherence factor 
measurement still provides a significant constraint 
in a combined measurement with two-body decays 
of the D. The precise effect depends on the exact 
value of the strong phase S^^; typically, without 
the K37r constraint from CLEO-c (but including 
the S^'^ constraint from CLEO-c), the precision for 
2fb-i data at LHCb from 0° ^ KK, tttt, Ktt, K37r 
modes is 10°, which improves to ^ 8° with the 
additional information from CLEO-c's coherence 
factor measurement in K37r [i^ . An equivalent study 
for Ktttt^ has not yet been performed, but the large 
value of Rktttt" suggests that we can expect a further 
significant improvement. 



6. Conclusion 



Charm physics has recently undergone a remark- 
able renaissance. Dalitz analyses in charm give access 
to magnitudes and phases of the charm decay ampli- 
tudes. This can be used to analyse light-meson reso- 
nance, where an improved description of broad S-wave 
resonances is of particular interest. Dalitz plots analy- 
ses in charm also play a significant role in charm mix- 
ing measurements, where Dalitz analyses provide the 
best constraints on x and important phase informa- 
tion that allows the translation of the x' , y' measure- 
ments in D — !■ Ktt into the unprimed parameters. And 
finally, in the decay chain — > D'^K*, and equiva- 
lently B° D°K*, phase information from the B de- 
cay is encoded in the D Dalitz plot, allowing a theoret- 
ically clean measurement of 7. This method provides 
currently the best constraint on 7. Input from the 
analysis of quantum-correlated DD pairs at CLEO- 
c will significantly reduce the model uncertainty in 
this group of measurements, in the D° — s- Kstttt mode 
as well as other 2, 3 and 4-body decay modes of the 
D°. This is of particular importance for future facili- 
ties such as the proposed Super Flavour Factory, and 
LHCb, which is due to start data taking in 2009, and 
where these uncertainties would soon be the limiting 
factor in the precision on 7. 
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